Bose-Einstein condensation and symmetry breaking 
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Adding a gauge symmetry breaking field — v\/V(aa + do) to the Hamiltonian of some simplified 
models of an interacting Bose gas we compute the condensate density and the symmetry break- 
ing order parameter in the limit of infinite volume and prove Bogoliubov's asymptotic hypothesis 
limy^oo(ao}/v / V : = sgnv limy_ 00 \J (a^ao) /V where the averages are taken in the ground state or 
in thermal equilibrium states. Letting v tend to zero in this equation we obtain that Bose-Einstein 
condensation occurs if and only if the gauge symmetry is spontaneously broken. The simplification 
consists in dropping the off-diagonal terms in the momentum representation of the pair interaction. 
The models include the mean field and the imperfect (Huang- Yang-Luttinger) Bose gas. An impli- 
cation of the result is that the compressibility sum rule cannot hold true in the ground state of the 
one-dimensional mean-field Bose gas. Our method is based on a resolution of the Hamiltonian into 
a family of single-mode (k = 0) Hamiltonians and on the analysis of the associated microcanonical 
ensembles. 

PACS: 03.75.Hh, 05.30.Jp 



I. INTRODUCTION AND RESULTS 

The relation between Bose-Einstein condensation 
(BEC) and a spontaneous breakdown of the gauge sym- 
metry of the system (shortly, symmetry breaking) is not 
entirely understood. Symmetry breaking is absent in 
textbook discussions of BEC of the ideal Bose gas and 
does not appear in the general definition of BEC by Pen- 
rose and Onsager [1] either. It is, on the other hand, 
a central notion of Bogoliubov's theory of superfluidity 
and superconductivity, dealt with extensively under the 
name of quasi-averaging [2] . Bose-Einstein condensation 
is the accumulation of a macroscopic number of particles 
in a single-particle state. In homogenous (translation in- 
variant) situation this state is <£>o( r ) = e^/yV where 
V is the volume of the cube on which the system is de- 
fined with periodic boundary conditions and is a real 
constant. The macroscopic occupation of (po means that 
the thermal (or ground state) average of its occupation 
number operator N increases proportionally to V as V 
tends to infinity, which can be shown to be equivalent 
to an off-diagonal long-range order (ODLRO). Hamilto- 
nians conserving particle number possess a continuous 
gauge symmetry related to the phase </>. For pair inter- 
actions a general Hamiltonian reads 



the periodic boundary condition. Averaging any o,k with 
eigenstates or the density matrix generated by H yields 
zero. (We disregard a possible accidental degeneracy of 
eigenstates belonging to different particle numbers.) For 
any real <fi changing all Ok simultaneously into e^^at 
(and a k into e l ^a k ) leaves the commutation relations 
and the Hamiltonian invariant. This gauge symmetry 
can be broken by adding a term which breaks particle 
number conservation, e.g. for real v and <f> the operators 
H — v\/V (e~ l( ^ao + e^ag) are no more gauge invariant. 
Still they are unitary equivalent for different </>, there- 
fore = can be taken without restricting generality. 
Now the thermal or ground state average (do) is non- 
vanishing, and we speak about a spontaneous breakdown 
of the gauge symmetry if the 'quasi-average' 
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(1.1) 



where e(k) = h 2 k 2 /2m, zx and u(q) = i>(— q) are real 
numbers, A k = a k ak, N = ^k, a k creates a boson in 
the single-particle state ^k(r) = (1/y/V) expjikr} and 
q, k, k' are allowed wave vectors, i.e. compatible with 



lim lim — =(an 

When this occurs, we are dealing with a continuum of in- 
finite volume pure states, each characterized by a phase 
</>, and the state we obtain by fixing v = in finite vol- 
umes and taking the limit V — > oo is the uniform mixture 
of pure states with < cj) < 2ir. 

By now the physical reality of symmetry breaking - 
at least in a non-translation-invariant situation - has be- 
come an experimental fact, manifesting itself in the in- 
terference of two condensates of different phases [3]. It 
is all the more embarrassing that there seems to exist no 
theoretical proof of its occurrence. There is Hohenberg's 
theorem on the absence of symmetry breaking in one and 
two dimensional Bose gases at positive temperatures [4- 
6]. In some other cases one can prove BEC by showing 
either the existence of ODLRO or the macroscopic occu- 
pation of a one-particle state. An example for the first 
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is the hard-core Bose lattice gas at half- filling [7] , for the 
second the trapped Bose gas [8]. However, the Schwarz 
inequality 



the order parameter 



|(ao)| < #o> 



(1.2) 



indicates that the absence of symmetry breaking does not 
automatically imply the absence of BEC, and vice versa, 
BEC does not obviously imply symmetry breaking. A 
closely related question is the validity of Bogoliubov's ap- 
proximation [9] which consists in replacing ao and by 
complex numbers c resp. c*. A partial answer was given 
by Ginibre a long time ago [10]. Ginibre considered this 
replacement in the Hamiltonian or in the grand-canonical 
density matrix, with or without the symmetry breaking 
field, and proved that in all cases by choosing c so that 
it maximizes the approximate pressure in finite volumes 
the result converges to the exact pressure in the limit of 
infinite volume. This suggests that symmetry breaking 
indeed takes place when there is BEC, and the inequality 
(1.2) saturates asymptotically. 

In this paper we will prove that for any v ^ 0, 



lim (a )/W = sgni/ lim tJ(N )/V 



(1.3) 



in the ground state and thermal equilibrium states of 
the mean-field Bose gas and of other models whose in- 
teraction is diagonal in momentum representation. The 
corresponding Hamiltonians are obtained from (1.1) by 
dropping the off-diagonal terms: 

i/ FD ^^e(k)iV fc - M 7V+^(iV 2 -iV) 

+ ^T7 E w ( k ' - V N * N * • ( L4 ) 

Here Vq = v(0) > 0, and for the sake of simplicity we 
consider only v(k) > 0. In the mean-field model, H MF , 
v(k) = for k 7^ 0. The imperfect Bose gas is defined 
by keeping the diagonal part of the (5 interaction, so that 
f(k) = vq for all k [11]. Note that only the interac- 
tion of the mean-field Bose gas is diagonal in coordinate 
representation. The other models may present nonin- 
tuitive features like anomalies in the spectrum of ele- 
mentary excitations, see below, and the Thouless effect, 
a jump-discontinuity in the condensate density as the 
chemical potential increases [12, 13]. The thermal equi- 
librium properties of the full diagonal model (1.4) were 
rigorously studied by Dorlas et al. [13]. 

The ground state of each Hamiltonian of the family 
(1.4) in the Fock space is \N = N — my) where my 
is the nonnegative integer closest to nV/v + \. In the 
limit of infinite volume the ground state energy density 
is — (i; /2) lim(my/V r ) 2 . We will study the Hamiltonian 

Hl D = H ¥D -vy/V[a + c^) (1.5) 
and prove that for v ^ the ground state energy density 
e = lim Eo/V = min{v n s 2 /2 — [is — 2\v\y/s} 

V^oo s>0 

= vqsI/2 - fiso - 2\u\^/so , (1.6) 



a = lim — =(a ) = sgnu 



and the condensate density 



Po = Jim t-(N ) 

V~*oo V 
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(1.7) 



(1.8) 



In (1.7) and (1.8) the averages are taken in the ground 
state of -ff FD . The minimizer sq coincides with the full 
particle density. It is positive for any value of fi, but 
goes to zero with a vanishing v if fj, < 0. The extension 
of Eq. (1.3) to positive temperatures will be obtained by 
proving the analog of Eqs. (1.6)-(1.8) for the ground state 
of certain microcanonical ensembles and by referring to 
the equivalence of micro- and grand-canonical ensembles. 
Equation (1.3) is not a statement about the existence of 
BEC. It holds even if the chemical potential is below 
its critical value when both sides vanish as v tends to 
zero. It however guaranties that BEC and a spontaneous 
breakdown of the gauge symmetry occur simultaneously. 

The result for the ground state is independent of the 
dimension, although iJ FD depends on it. In particular, 
there is symmetry breaking along with the Bose conden- 
sation in the ground state of the one-dimensional mean- 
field Bose gas. Because of the finite compressibility of 
this system, our finding is in conflict with an earlier ar- 
gument according to which there should be no symme- 
try breaking in the ground state of one-dimensional in- 
teracting Bose gases [14]. To track down the source of 
this contradiction and also by a general interest we con- 
sider here some implications of the results listed above 
on the properties of the excitation spectrum. The non- 
relativistic version of the Goldstone theorem states that a 
spontaneously broken continuous symmetry implies the 
existence of gapless excitations, i.e. eigenstates whose 
energy tends to zero with a vanishing momentum. Wag- 
ner's analysis [15] based on the Bogoliubov inequality and 
on Bogoliubov's 1 / g 2 -theorem concluded that the nature 
of the long- wavelength excitations was uncorrelated with 
symmetry breaking and was governed by the range of the 
interaction. Our conclusion for this particular family of 
interactions is similar. As it is shown below, there are 
gapless excitations in the mean-field model whose inter- 
action is of infinite range but there is a gap to excitations 
with a non-vanishing momentum in the spectrum of the 
other models. This seems to lead to another contradic- 
tion. In the ground state version of the 1 / g 2 -theorem, 



W > 7 



1 (^q-^IWI 2 



4 e(q) (N) + \v(ao)\\/V 



(1.9) 



Eq is the ground state energy and E^ is the energy of the 
lowest lying eigenstate of momentum fiq of the symme- 
try breaking Hamiltonian [16]. In the actual case (1.9) 
can be made stronger. The Hamiltonian i? FD commutes 
separately with each N^, k / 0. We shall prove that 
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for v 7^ the ground state is still in the invariant sub- 
space {Nk = 0|k 7^ 0}. Then the ground state average 
(A q ) — and from the derivation of Eq. (1.9) one sees 
that Eq can be taken to be the smallest eigenvalue among 
the eigenstates with A q = 1 and Ak = if k 7^ or q, 
i.e. a single particle carrying the total momentum, which 
can be much greater than the minimum eigenvalue among 
eigenstates of momentum fiq. In general, summing (1.9) 
over allowed q such that < |q| < qo, dividing by V 
and letting V tend to infinity while keeping q fixed, we 
arrive at 



P-Pa > 



1 



2(27r)° 



|q|<<jo 



1 | a | 2 lim v ^ 00 (£ q - E ) 
2 



e(q)p+ \ua \ 



1J dq 

(1.10) 



where p is the total particle density and d is the dimen- 
sion [17]. In the actual case p = po = \oto\ 2 . In Section 
IV we will derive the gap to one-particle excitations, 

E^-Eo=s^+^ + \ V \^±M 



+0(v 2 ) + 0(l/V), (1.11) 

valid for p > 0. The finite-size correction <~ 1/V is uni- 
form in v so that the limits V — ► oo and v — > are 
interchangeable. We can use this result in (1.10) and see 
what we obtain when v goes to zero. For the mean-field 
model hm„^ hmy^oo (-2^ — Eq) = £ (q)i therefore (1.10) 
holds true. For the other models the gap implies a di- 
verging integral in one and two dimensions and, thus, a 
violation of the inequality (1.10). The cause of the con- 
tradiction is our misuse of (1.9), which cannot be applied 
to if v(k) 7^ for some k ^ 0: These interactions 

are non-diagonal (nonlocal) in coordinate representation 
and therefore do not commute with /?k, the Fourier trans- 
form of the density operator, which is diagonal. This in- 
validates the /-sum rule [18], used in the derivation of 
Eq. (1.9). Turning to the controversy in connection with 
the mean-field model, we recall an improved version of 
(1.9) due to Pitaevskii and Stringari [14], 



(JVq>> 



1 |(a )| 2 1 1 |<a >| 2 



4 (p_qpq) 2 4 (A)S(q) 



where 5(q) is the static structure function. This formula 
is valid for all the models discussed in this paper and 
leads to 



/ ( 



2S(q) 



1 dq < . 



(1.13) 



We can conclude that q d / 5(q) — > as q — > in all models 
considered here. In the general case, when the /-sum rule 
holds true, (1.9) can be obtained also from (1.12). When 
both the /- and the compressibility sum rules [18] hold, 



lim q/S(q) > ^4m/h 2 Kp, 



(1.14) 



where k is the ground state compressibility. In the mean- 
field Bose gas k = l/vop 2 < oo. Thus, (1.14) cannot hold 
true in one dimension and because the /-sum rule is valid, 
the compressibility sum rule must fail in the ground state 
of the one-dimensional mean-field Bose gas. 

Note that the positive temperature variant of inequal- 
ity (1.9), with 2k B T replacing Eq - E , cf. [2, 16], also 
relies on the commutability of /?k with the interaction. 
Thus, it cannot be applied to prove the absence of sym- 
metry breaking for iJ FD in one and two dimensions at 
T > if v(k) is nonvanishing for nonzero k. Rightly so, 
because in these systems, e.g. in the imperfect Bose gas 
[11], there can be BEC at T > even in one dimension, 
see also [12, 13]. 

All our results arc obtained by representing the Hamil- 
tonian as a Jacobi matrix in microcanonical ensembles or, 
equivalently, by studying some tight-binding Schrodinger 
equations on the spectrum of A . In Section II we show 
that Ak = for all k 7^ in finite volume ground states. 
Formulas (1.6)-(1.8) will be derived in Section III. The 
analog of (1.6) for the ground state of microcanonical en- 
sembles is presented there, that of (1.7) and (1.8) in Sec- 
tion IV which also contains the derivation of Eq. (1.11). 
In Section V we discuss the extension of the results to 
positive temperatures. The paper ends with a Summary. 



II. GROUND STATE IN FINITE VOLUMES 

The Hamiltonians commute with each Ak, k ^ 0. One 
can therefore decompose them into operators acting on 
invariant subspaces with fixed N^. Let n = {nkjk^o 
denote terminating sequences of nonnegative integers, i.e. 
Il n ll = Sk#o nk < 00 ■ Trien 



H, 



FD 



)n#(n) 



(2.1) 



where H (n) is the restriction of onto the subspace 
{Ak = ^k|k 7^ 0}. Because v > 0, H™ has a discrete 
point spectrum with +oo as the only accumulation point, 
Tr exp{— (3H™} < oo for (3 > 0, and the same holds 
for each H(n). In this section we prove the following 
assertion. 



(1-12) In any finite volume the ground state of 



FD 



subspace {Ak 
H(0). 



i the 

0|k 7^ 0}, i.e. it is the ground state of 



This property is specific to the diagonal models and 
v(k) > 0. A general interaction does not commute with 
Ak. This knowledge about the ground state can but need 
not be used for the derivation of the formulas (1.6), (1.7) 
and (1.8). It plays, however, a role in the interpretation 
and computation of Eq — E in Eq. (1.11). 

One has to prove that for each n 7^ the ground 
state energies of H(n) and H(0) satisfy the inequal- 
ity E [H(n)] > £b[-H"(0)]. This is immediately seen if 
||n|| > 1 + 2pV/v . The operator H(n) - H(0) is di- 
agonal and in this case it has positive elements, which 
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implies the result. In the case ||n|| < 1 + 2/j.V/vq (which 
can occur only if /j, > 0), however, it may have both 
positive and negative elements and we need a more elab- 
orate argument. We shall interpose between -ff(n) and 
H(0) three auxiliary operators, if(||n||), X(||n||,0) and 
K (||n||, 1) and show the following order of ground state 
energies: 

E [H(n)} > E [K(\\n\\)} = E [K(\\n\\, 0)] 

>E o [K(\\n\\,l)]>E [H(0)}. (2.2) 

Define 

Kip) = ^(7Vo+p) 2 -( M + ^) (N +p)-vVV(a +a* ) 

(2.3) 

where p is a real number. In particular, K(0) = H(0). 
The operator inequality 



H(n)>K(\\n\\ 



(2.4) 



holds true for n^O, because H(n) — K(\\n\\) is a diagonal 
operator with entries > e(k)rik > 0. For the ground 
state energies (2.4) implies 



Eo[H(n)]>Eo[K(\\n\\ 



(2.5) 



Let p be a nonnegativc integer. Using the basis 
{\N — m)}m = o °f -ZVo-eigenstates, K{p) can be rep- 
resented by a semi-infinite tridiagonal matrix (Jacobi 
matrix), denoted also by K(j>). The diagonal matrix 



elements satisfy K(p) mm = K(0) m+Ptm+p for m > 
0. The off-diagonal elements are the same for all p, 
K{p)m-i,m = K(p)m, m -i = -v^fmV if to > 1, and 
otherwise. The matrices K(p) act in the Hilbcrt 

space &(N) = {(VV))o°jEjV'(m)| 2 < °°}- We ex- 
tend them without changing the notation onto t^ilj) = 
{(V>(™))-JE m |^(to)| 2 < ex.} by setting K{p) mn = 
if cither m or n is negative. Then t%(N) and ^(2\N) are 
invariant subspaces of each K(p) which acts on the first 
space as the original operator and on the second as the 
null operator. If Eo[K(p)] < 0, the ground state energy 
of the extended operator is the same as that of the orig- 
inal one. Let U p be the left shift by p, i.e. for ip € ^2(^)5 
(U p tp)(m) = tp(m+p). Then 



{U^K(p)U p ) m 


K(0) r 



— K{p)rn—p,m—p 
, TO < p — 1 

n , m>p 



(2.6) 



and 



(U-iKWUpU-!^ 



-VyJ (to — p)V , m>p+l 



K{p)m— 1— p,m— p 
, TO < p 



(2.7) 



Define K(p,0) = U p 1 K(p)U p . The operator U p is uni- 
tary, thus E [K(p, 0)] = E [K(p)]. 

For the definition of K (p, 1) we split the matrices 
K(p 7 0) and K(0) into diagonal D and off-diagonal Q 



parts, K(p,0) = D{p) + Q(p) and K(0) = D(0) + Q(0), 
and set K (p, 1) = D(p) + Q(0). Furthermore, we intro- 
duce an interpolating matrix function 

K(p, t) = tK(p, l)+(l-t)K(p, 0) = D(p)+t[Q(0)-Q( P )} . 

(2.8) 

By the Hellmann-Feynman theorem 

9 -E [K(p,t)} = f> t , [Q(0) - Q( P Mt) 



dt 



Mm)[Q(0) - Q(p)} m nMn) (2.9) 



m.n>0 



where ipt is the normalized ground state of K(p,t). 
Choose v to be positive (changing the sign of v amounts 
to a unitary transformation with a diagonal matrix of 
elements (-!)"')• Then 



[0(0) - Q{p)]m-l,m = [0(0) - Q{p)} 



m.m—1 



< 



for to > 1 implying (through the variational principle 
or the Perron- Frobenius theorem) Vt( m ) > for to > 
0. Thus the derivative (2.9) is negative, which proves 
E o [K(p,0)}> E [K(p,l)}. 

At last, we note that K(p, 1)-K(0) = D(p)-D(0) > 
if p = ||n|| < 1 + 2[iV/v because 



[D(p)-D(0)] r 



K(0) r 



> 







< TO < p — 1 

otherwise . 

(2.10) 

Therefore £ [K(||n||, 1)] > E o [K{0)} indeed. Here 
i?o[i^(0)] is the ground state energy of the extended 
operator. If E [H(0)] was positive, we would obtain 
E o [K(0)} = 0. However, E o [H(0)} < E [H FU ] < can 
be seen by using \N = No = my) as a variational wave 
function, so that E [K(0)] = E [H(0)}. 



III. GROUND STATE IN INFINITE VOLUME 

Our aim in this section is to derive Eqs. (1.6)-(1.8). 
Incidentally, we shall obtain the ground state energy den- 
sity for H(n), cf. (2.1). 

Each H (n) is an unbounded operator which is bounded 
below. A common lower bound is that of H™ , 



H, 



FD 



> H 



> - 



FD 



\v\(V + N ) 



T 



(J, 



v 



+ 



2V 



+ W 



V. (3.1) 



Here we have made use of ±W(a + al) <V + N a which 
comes from (W ± ag)(-\/U ± ao) > 0. 

To obtain a good variational upper bound on 
E [H(n)], note that i? FD of Eq. (1.4) has a decom- 
position analogous to (2.1), H FB = n ffdiag(n) where 
iJ(iiag( n ) is the diagonal part of H(n). Let |iV = J(n)) 
be the ground state of £?diag(n) (J(0) = my). Then with 



t/j = V 



- 1/4 E i*o 

m— J(n) 



m) 
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as a variational wave function we obtain 

E [H(n)} < E [H diag (n)} - 2\v\^J(njv + 0(VV) . 

(3.2) 

This upper bound will turn out to be the precise asymp- 
totic result up to order \v\. 
Let us introduce 



(n) = -^^ nk Vv(n) = ^^2v(k)n k 



V 



k^O 



2V 2 



n k n k / 



(3.3) 



All three are real nonnegativc numbers, vanishing if 
II nil = 0. With them 



H(n)/V = h v {x Vl y Vl z v ) 



( , v ° \ ( Na 



V 



+ x v 



+W— - ^( fl o + «o) + z v ■ 

(3.4) 

As in the previous section, we consider the matrix of 
each hv in the basis of the A^o-eigenstates and extend 
it to £2(2) with zero elements for negative indices. The 
nonvanishing matrix elements of hv are 

hv{xv,Vv,Zv)jj = 

?(M 2 -(<^)(M+»4+- 

(3.5) 

and 



h v (x v ,yv,zv)j,j+i = h v (x v ,yv,zv)j+i,j = -v 



3+1 
V 

(3.6) 



with j > 0. 

Each hv is still an unbounded operator. We can zoom 
in on different parts of its spectrum by employing dif- 
ferent Independent right shifts in £2^) and taking the 
limit V — > 00. (Left shifts would only generate the zero 
matrix.) Suppose that xv — > x, yv — > V and zy — > z. 
The limits x, y and z may not be arbitrary nonnega- 
tivc numbers, e.g. y = for the mean-field model and 
y = vqx for the imperfect Bose gas. Choose a sequence 
Sv of positive integers such that sv — Sv/V — > s. For 
the matrix elements the shift corresponds to a change of 
variables, j = m + Sv, yielding hv{xv, yv, Zv,sv)mm' — 
hv(xv,yv,Zv)m+S v ,m'+Sv When V tends to infinity 
the limit of hy(xv , yv , Zy, sy) mm i exists and defines a 
doubly infinite Jacobi matrix h(x, y, z, s) with nonvan- 
ishing elements [19] 

h(x, y, z, s) mm = x(x, y, z, s) 

= y( s + a; ) 2 -v(s + x) +ys + z (3.7) 



and 

h(x, y, z, s) m , m+ i = /i(a;, y, z, s) m+ i, m = -v\fs . (3.8) 

Formally h(x,y, z, s) is a one-dimensional tight-binding 
Schrodinger operator with a constant potential (3.7), 
governing the motion of a particle of mass ~ l/\i/\^/s 
along the spectrum of Nq. Physically h(x, y, z, s) is asso- 
ciated with the infinite volume limit of a microcanonical 
ensemble of bosons characterized by four intensive vari- 
ables: the condensate density s, the density x and energy 
density —\xx + vqx 2 /2 + z of the uncondensed particles 
and the density of the interaction energy, (y + vqx)s, 
between the condensate and the uncondensed particles. 
Because of the symmetry breaking field, the energy den- 
sity of the condensate is still an operator of diagonal and 
off-diagonal matrix elements —/is + v a s 2 /2 and —v^fs, 
respectively. 

For s > the spectrum of h(x, y, z, s) is purely abso- 
lutely continuous. As a set, it is the interval 

[x(x, y, z, s) - 2\v\y/s, x(x, y, z, s) + 2\v\sfs] 

that can be found by solving the generalized eigenvalue 
equation hip = etp with the usual exponential ansatz 
ip(m) = cxp(ima) and choosing any real a in the interval 
[0,7r]. The spectral point corresponding to a is 

e(x, y, z, s, a) = x( x , y, z, s) — 2v\fs cos a . (3.9) 

In particular, the lower edge of the spectrum is 

9xyz(f*, s) = ys 2 - (fi - v x - y)s - 2\v\y/s 

-flX+ yV +Z (3.10) 

with two linearly independent solutions ip(m) = 1 and 
ip(m) = m. The ground state energy density of H(n) is 
obtained by minimizing the lower edge value with respect 
to s, 



lim E [H{n)]/V = < 



rnm.g xyz (iJL,v,s) 

s>0 



xyz ^) 
= 9xyz{H,V,S X y(ll, v )), (3- 11 ) 



where the minimizer s xy does not depend on z [20]. In- 
specting the expression (3.10) we see that whenever v ^ 
0, s xy > for any value of [i. Hence, dg xyz ()i, v, s)/ds = 



at s — s xy , that is, 



is a positive root of a cubic 



polynomial. 

We have seen in the previous section that in the overall 
ground state n = 0, therefore x = y = z = provides 
the overall ground state energy density: 



eo([A, v) — min{w s 2 /2 — /j,s — 2\v\y/s} 

s>0 

= mm g (fj,,v,s) = g (fi,v,s (tJ,,v)) . 

S>0 



(3.12) 



Equation (3.12) agrees with Eq. (1.6). We could have 
found eo without knowing the ground state in finite vol- 
umes, by minimizing e xyz over nonnegative values of x, 
y and z. 
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To obtain Eq. 
tions: 



(1.7) we make the following observa- (i) Choosing v ^ 0, 



dv 



dgoij^u, s) 
dv 



+ 



<9sqQ, v) 
dv 



dg (n,v, s) 



ds 



dg (n,v,s) 

dv 



= -2sgniV^- (3.13) 



s=s 



(ii) The ground state energy Eq (/x, v) of H„ is a con- 
cave (even) function of v. Concavity comes from that of 
the thermodynamic potential —fl~ 1 ln Tr exp{— (3H™} 
by letting the temperature go to zero. (E (fi, v) is also 
an analytic function of v in finite volumes; the singularity 
of Eq/V at v = develops only as V tends to infinity.) 

(iii) The Hcllmann-Feynman theorem yields the ground 
state averages 



(ao> = (a* ) 



1 dE Q (/i,v) 
2s/V dv 



(3.14) 



(iv) According to an inequality known as Griffiths' lemma 
[21] and applied here to the sequence of concave functions 
E (fi, v)/V converging to e (n, v), 

de (n,v + 0) <liminf 1 dE (ii,v) 



dv 



V^oo V 



dv 



\ dE Q {n,v) de (n,v-0) 
<hmsup — < (3.15 

y^oo V dv dv 

implying equalities in points where e is differentiablc. 
With (3.14) and (3.15) for v ^ 



a 



lde (n, v) 
2 



dv 



(3.16) 



and substituting (3.13) into this equation we find Eq. 
(1.7). 

Eq(h, v) being a concave function also of /z, Eq. (1.8) 
can be proved from 



d[i 



dgo(p,v, s) 



d/j, 



= -so 



(3.17) 



by an identical argument. Because in the Hamiltonian 
/i is coupled to N and not to Nq, we can use the result 
of the former section, that the ground state expectation 
value of Nk vanishes if k ^ 0. Therefore 



Po 



de (p, v) 
d\i 



(3.18) 



which together with (3.17) yields Eq. (1.8). 

For it > by solving the cubic equation we find 



\ 3 r 



(3.19) 



with q = —fi/3v and r = \v\/2v . Expansion of ^/so up 
to first order in \v\ yields 



v 2/i 




(3.20) 



where /,(/x, v) are smooth functions away from zero and 
bounded at v = 0, and ji and are positive. 



IV. GAP TO ONE-PARTICLE EXCITATIONS 
WITH A NONVANISHING MOMENTUM 

Let us return to Eq. (3.11). The minimum of 
9x V z{n, 0, s) is attained at s = (fi — fi c (x, y))/v which is 
positive if \i > fi c (x, y) = y + v x. This is the condition 
to have a nonvanishing Bose condensate at v = in the 
ground state of H(n) with n chosen so that xy(n) — ► x 
and yv( n ) — > y, cf. Eq. (3.3). If v ^ and \x < fj, c (x,y), 
the ground state condensate density s xy is still positive 
but tends to zero together with v. Now the full density 
is x + s xy . Properties (i)-(iv) listed in the former section 
remain valid, thus the order parameter and the conden- 
sate density in the ground state of H(n) will respectively 
be 



and 



>-xy 



Pxy 



Ide 
2 



xyz 



dv 



sgn vJs, 



xy 



de xyz {[i,v) 
dfj, 



Jxy 



(4.1) 



(4.2) 



We conclude that the relation 



Pxy(p,v) 

holds true for any fx, provided that v =/= 0. For /i > 
fi c {x, y), s xy is given by the right-hand side of Eq. (3.19) 
with q = —(pi — [i c (x, y))/3vo and r = \v\/2vo- 

Substituting s xy into Eq. (3.10) we obtain the ground 
state energy density for H(n). To prove the formula 
(1.11) this general expression is not needed. We set 
n q = 1 and rik = for k ^ q. Then 

H(n) - H(0) = e(q) - M + [v + v(q)]N /V . (4.3) 
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Averaging with the ground states of H(n) and H(0), re- 
spectively, 

e(q)-J*+[ifo + v(q)]<JVb) q /V r 

< £ q - So = £o[#(n)] - E [H(0)} 
<£(q)-p+[v +v(q)](N ) /V 

(4.4) 

by the variational principle. In the present case x = y = 
z = 0, so both (N )q/V and (N ) /V equal p + O(l/V). 
For p > we insert po from (3.20) and obtain Eq. (1.11). 



V. EXTENSION TO POSITIVE 
TEMPERATURES 

Our analysis of microcanonical ensembles, especially 
equations (4.1) and (4.2) make it possible to extend the 
result (1.3) to thermal equilibrium states. Below we out- 
line the argument leading to this conclusion. 

In the microcanonical ensemble {x, y, z, s} the conden- 
sate density is s, while the order parameter operator 
limy^oo ciq/VV is represented by a matrix o with ele- 
ments o m m+ i = ^/s, — oo < m < oo and otherwise. 
We can compute the mean value of o in a generalized 
eigenstate tp of h(x, y, z, s) by the formula 



(o)v, = lim 

M^oo 



M<m<M 



ip*(m)(atp)(m) 



-M<m<M 



(5.1) 



Consider ip ±a (m) — e ±tma , Va,c( m ) = cosma and 
ipa,s(m) = sin ma, all belonging to the spectral point 
(3.9). Now (o) v , ±0 . = e d ' 



" ±ta y/s, while 



(°)V-a,c = ( )l/Vs = COS OLtJ~S. 



(5.2) 



We want to treat the above mean values as limits of 
averages taken with eigenvectors of . By a gen- 
eral property of one-dimensional second order difference 
equations, all eigenvalues of H(n) are nondegenerate and 
therefore the eigenvectors are real. This property is in- 
herited by H™ if we disregard any possible accidental 
coincidence of eigenvalues belonging to different n. So 
considering separate eigenvectors the relevant relation is 
(5.2) in the sense that for any sequence of eigenvectors 
ipi, the volume increasing with I, for all possible limits 



lim 



1 (i/ji\ao\ipi} 



1 / MjjVoM) 
cos a lim — = — 

W V lM>l>i) 



1/2 



(5-3) 

where a can assume any value. Note also that computing 
the average with the help of a density matrix we have to 
take a trace and then find the same outcome (5.2) or 
(5.3) whether we use the complex basis tp± a or the real 

one, ^ a ,c,1pa,s- 

Unless a = or w, Eq. (5.3) is in conflict with Bogoli- 
ubov's hypothesis, supposing an asymptotic equality in 
the Schwarz inequality (1.2). The fact that Bogoliubov's 



suggestion proves nevertheless to be correct and we ob- 
tain Eq. (1.3) in thermal equilibrium is due to the strong 
equivalence of the grand-canonical and a suitably chosen 
microcanical ensemble: fixing /3 = l/fc^T, p and v and 
taking the thermodynamic limit, the distribution of the 
intensive variables, e, x, y, z and s among them, becomes 
a Dirac delta. Let the energy density be ep^ v . Although 
there is a continuum of the parameter sets {x,y,z,s} 
yielding e(x, y, z, s, a) = ep^ v for different values of a, 
cf. (3.9), because of the degenerate limit distribution 
of these variables only a unique set will asymptotically 
dominate. With this set ep^ u is attained as the mini- 
mum of e(x, y, z, s, a) with respect to a and s because 
the density of states of each H(n) has a maximum at the 
bottom of the spectrum of H(n). Minimization singles 
out a — or ir, depending on the sign of v; see Eq. (3.10). 
In conclusion, ep^ = e xyz (p, v) and the thermodynamic 
limits of (a Q )p l _ lu /W and (N }p f _ l „/V are a xy resp. p xy , 
cf. Eqs. (3.11), (4.1) and (4.2), with x = x((3,p,if), 
y = y(p,p,i>) and z = z(f3,p,v). This proves the va- 
lidity of Eq. (1.3). Note that p c ((3), the critical chemical 
potential of the grand-canonical ensemble satisfies the re- 
lation 

HciP) = p c (x(/3, /i c (/3), 0), y(/3, p c (/3), 0)) 

= y((3, p c ((3), 0) + vox(0, p c ((3), 0). (5.4) 

Recently Pule and Zagrebnov [22] computed the ther- 
modynamic limit of the pressure of the mean-field Bosc 
gas in the presence of a gauge-symmetry breaking field. 
Their formula reads 

p(p, v) = Jim In Tr e _/3H " F = ^-p 2 (p, v) 



v^oo (3V 

+Pideal(M- Vop{p.,v)) + 



v p(p, v) - p 

where p(p, v) is the unique solution of the equation 

,.2 

P = Pidcalip - V 0P) + 



(5.5) 



(v Q p-p) 2 ' 



(5.6) 



Here Pideai(p) and Pidea,i(p) are respectively the pressure 
and density of the ideal Bose gas for v = 0. They are well- 
defined if the chemical potential is negative. The authors 
prove that for v ^ 0, p — v p{p, v) < indeed. Because 
the pressure in any finite volume is a convex function 
of i>, we can use Griffiths' lemma to obtain the order 
parameter from the relation a xy = a x o = [l/2)dp/dv. 
We find a x0 — vj (v p(p, v) — p), thus 

Pxo = alo = v) - Pideai(M - v p(p, v)). (5.7) 

Since p(p, v) is the total density, we can identify the den- 
sity of the uncondensed particles as 

x = Pid ea \{p-v p{p,v)). 

If p < p c (/3) = foPideai(O), cx x o goes to zero with v. in this 
case p — vop(p, v) converges to a negative value, p(/x,0) 
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is the unique solution of the equation p = pn elx \{p — 
v p). This latter holds true at p = p c (/3) when p c ((3) — 
vop(p c {(3), v) tends to zero with a vanishing v, but slower 
than v. If p > p c (/3), the convergence to zero becomes 
linear and thus the limit of a x o will be nonzero. 

VI. SUMMARY 

In this paper we have rigorously established the equiv- 
alence of Bose-Einstein condensation and a spontaneous 
breakdown of the gauge symmetry in some simplified 
models of interacting bosons. In these models one retains 
only that part of the pair interaction diagonal in mo- 
mentum representation. Equivalence has been obtained 
in the form of Eq. (1.3), that we found to hold true for 
any nonzero value of the symmetry breaking field, not 
only in the limit of a vanishing field. As a by-product, 
we have disproved the compressibility sum rule in the 
ground state of the one-dimensional mean-field Bose gas. 
Our method was to study certain microcanonical ensem- 



bles with a fixed density and energy density of uncon- 
densed particles. In the thermodynamic limit these en- 
sembles are still described by a family of operators, so- 
called Jacobi matrices, known from the one-dimensional 
tight-binding electron theory. The limit of infinite vol- 
ume could be done on the many-body Hamiltonians be- 
cause of their commuting with the occupation number 
operators of nonzero momentum states; this made it pos- 
sible to reduce the problem to that of a single mode, 
k = 0. The result was first obtained in the ground state 
of the microcanonical ensembles and then extended to 
thermal equilibrium states by arguing with the strong 
equivalence of grand- and microcanonical ensembles. 
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